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•  A  TEG  model  is  proposed  and  validated  under  transient  and  steady-state  conditions. 
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Thermoelectric  generators  (TEGs)  make  use  of  the  Seebeck  effect  in  semiconductors  for  the  direct  conver¬ 
sion  of  heat  into  electrical  energy,  being  of  particular  interest  for  high  reliability  systems  or  for  waste  heat 
recovery.  Although  several  TEG  models  can  be  found  in  the  literature,  many  of  them  not  offer  a  theoret¬ 
ical  solution  because  they  are  based  on  steady-state  solutions  or  they  are  assuming  fixed  parameters  as 
boundary  conditions.  Consequently,  to  assess  and  optimize  thermoelectric  generators  in  real  applications 
a  numerical  transient  simulation  tool,  which  takes  into  account  the  whole  energy  system,  is  mandatory. 
For  that  purpose,  a  new  TRNSYS  type  is  developed.  This  TEG  component,  which  can  be  used  as  a  design 
tool,  is  presented  in  this  paper  and  validated  using  experimental  data. 

The  results  show  that  the  proposed  component  is  able  to  cope  with  both  thermal  and  electrical  dynam¬ 
ics.  The  comparison  between  theoretic  and  experimental  results  has  approved  the  reasonability  of  the 
new  component.  The  normalized  root  mean  square  errors  are  3.53%  and  2.33%  for  temperature  difference 
between  hot  and  cold  sides  and  electrical  output  power,  respectively. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

A  thermoelectric  power  generator  is  a  solid-state  device  that 
provides  direct  energy  conversion  from  thermal  energy,  due  to  a 
temperature  gradient,  into  electrical  energy  based  on  Seebeck 
effect.  Also,  they  can  work  in  reverse  and  use  electrical  energy  to 
create  a  temperature  gradient  for  cooling  or  heating  applications. 
The  absence  of  moving  parts,  wide  range  of  operating  tempera¬ 
tures,  scalability,  and  modular  capabilities  makes  thermoelectricity 
attractive  for  a  wide  variety  of  applications,  such  as  power  for 
remote  control  and  monitoring  of  oil  or  gas  pipelines  and 
production  facilities,  automotive  waste  heat  recovery,  power  for 
navigational  aids,  spacecraft  radioisotope  power  supply,  telecom¬ 
munications  systems  and  cathodic  protection,  and  other  energy 
recovery  processes  [1-6],  Thermoelectric  devices  have  relatively 
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low  efficiencies  but  there  have  been  recent  advances  in  thermo¬ 
electric  materials  potentially  opening  the  door  to  new  power 
applications  [7,8],  As  material  advancements  continue  and  a  wider 
range  of  power  generation  applications  will  be  considered,  module 
and  system  level  modelling  becomes  critical  for  the  design  of  the 
next  generation  of  thermoelectric  systems. 

Although  several  models  for  TEG  modules  can  be  found  in  the 
literature  [9-16],  many  of  them  not  offer  a  theoretical  solution 
because  their  governing  equations  are  based  on  steady-state 
solutions  or  they  are  assuming  fixed  temperatures  as  boundary 
conditions  at  both  sides  of  the  TEG.  Moreover,  almost  none  of  them 
have  studied  the  transient  effects  of  load  resistance.  Complete 
transient  analyses  are  seldom  presented  and  just  a  few  have 
already  provided  a  complete  mathematical  solution  of  the  heat 
conduction  equation  for  TEG  devices  [13-15], 

This  study  attempts  to  fill  the  existing  gap  in  the  simulation  of 
thermoelectric  generation  through  the  development  of  a  new 
component  that  can  be  used  in  TRNSYS  software.  TRNSYS  [17], 
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Nomenclature 

T 

temperature  (K) 

Abbreviation 

y 

axis  coordinate  (m) 

TEG 

Symbols 

thermoelectric  generator 

Subscript 

time  (s) 

n 

number  of  thermocouples 

hp 

hot  plate 

a 

Seebeck  coefficient  (V/K) 

tCl,2 

thermal  compound 

a 

electrical  resistivity  (C 2  m) 

ce 

ceramic  substrate 

e 

length  (m) 

cp 

cold  plate 

A 

cross  sectional  area  (m2) 

cb 

water-cooling  block 

P 

material  density  (kg(m3)) 

h 

hot  side 

S 

specific  heat  capacity  (J/kg  K) 

c 

cold  side 

X 

thermal  conductivity  (W/m  K) 

P 

p-type  semiconductor 

K 

thermal  conductance  (W/K) 

n 

n-type  semiconductor 

C 

thermal  capacity  (J/K)) 

in 

input 

h 

convective  heat  transfer  coefficient  (W/m2  K) 

L 

load  resistance 

Q 

heat  rate  (W) 

Coo 

cooling  water 

V 

electric  voltage  (V) 

oc 

open  circuit 

I 

electric  current  (A) 

max 

maximum 

P 

n 

R 

electric  power  (W) 
efficiency  (%) 
electric  resistance  (D) 

ss 

steady-state 

developed  at  the  University  of  Wisconsin,  is  a  transient  systems 
simulation  program  with  a  modular  structure.  It  recognizes  a  sys¬ 
tem  description  language  in  which  the  user  specifies  the  compo¬ 
nents  that  constitute  the  system  and  the  manner  in  which  they 
are  connected.  The  TRNSYS  library  includes  many  of  the  compo¬ 
nents  commonly  found  in  thermal  and  electrical  energy  systems, 
as  well  as  component  routines  to  handle  input  of  weather  data 
or  other  time-dependent  forcing  functions  and  output  of  simula¬ 
tion  results.  The  modular  nature  of  TRNSYS  gives  the  program  flex¬ 
ibility,  and  facilitates  the  addition  to  the  program  of  mathematical 
models  not  included  in  the  standard  TRNSYS  library.  TRNSYS  is 
well  suited  to  detailed  analyses  of  any  system  whose  behaviour 
is  dependent  on  the  passage  of  time.  It  has  become  reference 
software  for  researchers  and  engineers  around  the  world.  Main 
applications  include:  solar  systems  (solar  thermal  and  photovoltaic 
systems),  low  energy  buildings,  HVAC  systems,  renewable  energy 
systems,  cogeneration,  fuel  cells  and  more. 

Full  integration  of  the  new  component  within  the  TRNSYS  sim¬ 
ulation  package  is  another  advantage  of  this  study,  which  makes  it 
more  applicable  for  designers  in  both  the  design  and  commission¬ 
ing  of  waste  energy  harvesting  systems.  The  incorporation  of  the 
TEG  model  into  standard  TRNSYS  library  will  allow  users  to  simu¬ 
late  TEGs  coupled  to  many  different  thermal  facilities  and  applica¬ 
tions  (e.g.  solar  thermal  systems,  geothermal  systems  and  more). 
Because  of  the  mentioned  benefits,  TRNSYS  is  chosen  to  implement 
the  new  component. 

The  objective  of  this  work  is  to  develop  a  computational  model 
capable  to  simulate  the  electro-thermal  dynamics  of  a  TEG  system 
into  TRNSYS  software.  The  new  component  is  fully  analysed  and 
validated  under  steady  and  transient  states  with  data  obtained 
from  the  experimental  setup. 


2.  TEG  modelling 

A  TEG  device  in  essence  is  a  thermopile  composed  of  a  number 
of  p-  and  n-type  semiconductor  pairs  connected  electrically  in 
series.  Heat  carries  the  majority  carriers  from  one  junction  to  the 
other  producing  a  current  and  voltage.  By  placing  many  PN  couples 
in  series  electrically  and  in  parallel  thermally,  a  TEG  module 


generates  an  open-circuit  voltage  proportional  to  the  temperature 
differential  across  the  elements. 

According  to  Seebeck,  Peltier  and  Thomson  effects,  the  TEG 
power  generation  depends  mainly  on  both  temperature  differen¬ 
tial  across  the  semiconductor  elements  and  electrical  load 
connected  (i.e.  changes  of  the  temperature  difference  lead  to 
changes  of  the  open-circuit  voltage,  and  the  TEG  output  perfor¬ 
mance  will  change  as  well).  Therefore,  the  unsteady-state  heat 
transfer  model  is  built  to  solve  the  TEG  temperature  distribution. 
Thereafter,  the  TEG  power  and  performance  outputs  can  be 
determined. 

Similar  with  Gou  [12,21],  Kim  [20],  Nguyen  [14]  and  Rodriguez 
[19],  the  basic  configuration  of  the  thermoelectric  study  follows 
the  general  design  illustrated  in  Fig.  1 . 

A  thermoelectric  module  is  sandwiched  between  electrical 
heater  and  water-cooled  heat  exchanger  that  can  be  treated  as  a 
convective  heat  transfer  system.  The  TEG  consists  of  certain  num¬ 
ber  of  semiconductors,  corresponding  solder  layers,  conducting 
strips  and  thermal  insulation  material.  In  order  to  electrically  iso¬ 
late  the  module,  a  ceramic  substrate  is  located  both  sides  of  the 
TEG.  A  thermal  compound  increases  the  thermal  conductivity  of 
the  interfaces  by  filling  microscopic  air-gaps  between  ceramic  sub¬ 
strates  and  aluminium  plates  and  TEG  and  ceramic  substrates, 
respectively.  At  the  same  time,  hot  and  cold  plates  allow  to  intro¬ 
duce  temperature  sensors  on  the  experimental  set  up. 


2.2.  Assumptions  and  boundary  conditions 

In  the  analytical  assessing,  the  thermal  conductivity,  electric 
conductivity,  and  Seebeck  coefficient  of  semiconductors  are  all 
assumed  to  be  temperature-independent,  and  adiabatic  boundary 
conditions  were  supposed  on  the  side  surfaces  of  TEG  element. 
Thermal  conductivities,  electrical  resistivity,  and  specific  heat 
capacities  of  non-thermoelectric  materials  are  supposed  constants 
within  the  operating  temperature  range. 

Due  to  the  fact  that  this  study  is  focused  on  low-temperature 
system,  the  Thomson  effect  could  be  neglected  [17],  Besides,  initial 
temperature  of  the  system  is  equal  to  ambience  temperature,  this 
means  that  temperatures  at  both  sides  of  the  TEG  element  are  the 
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Fig.  1.  One-dimensional  heat  transfer  model  of  thermoelectric  generator. 


same,  and  consequently  initial  output  data  is  zero.  Finally,  an 
instantaneous  response  to  electrical  load  transients  is  considered 
[13]. 

The  heat  leakage  through  the  solder  layer  and  conducting  strips 
are  neglected.  The  flow  of  heat  and  current  in  the  system  are 
assumed  one-dimensional.  All  materials  are  assumed  homoge¬ 
neous  and  isotropic.  Finally,  heat  losses  due  to  radiation  and  trans¬ 
verse  convection  are  neglected. 

2.2.  Governing  equations 

In  thermoelectric  devices,  both  refrigeration  and  generation 
applications,  the  internal  and  external  nodal  temperatures  of  TEG 
cannot  be  determined  analytically  without  knowing  the  heat  flux 
due  to  the  thermoelectric  effects.  The  reasons  are  the  nonlineari¬ 
ties  of  the  thermoelectricity  and  heat  transfer  equations. 

The  model  solves  the  nonlinearities  using  finite  difference  and 
Newton-Raphson  methods,  which  calculates  the  temperature  at 
different  nodes  separated  in  space  by  a  discrete  distance.  In  the 
transient  state,  the  temperatures  of  these  points  are  calculated  at 
discrete  periods  of  time  and  the  temperatures  for  all  the  nodes 
are  recalculated  at  the  end  of  this  time  interval.  Using  the  implicit 
finite  difference  method,  the  values  of  heat  flux  can  be  determined 
using  the  values  of  the  temperatures  of  the  time  step  before. 

According  to  heat  conduction  theory,  the  one-dimensional 
unsteady  heat  conduction  equation  of  non-thermoelectric  ele¬ 
ments  can  be  shown  as: 


PtSi 


m 

dt  ~  M  dyf 


(1) 


where  p,  S  and  X  are  density,  specific  heat  capacity  and  thermal 
conductivity,  respectively,  t  and  y  are  time  and  axis  coordinate, 
respectively.  The  subscript  i  determines  the  heat  transfer  related 
to  node  i.  As  found  in  the  referenced  literature  [13,19],  the  ther¬ 
mal-electrical  analogy  is  useful  in  the  analysis  of  complicated 
unsteady  heat  transfer  problems,  which  can  be  understood  by  cre¬ 
ating  an  electric  circuit  like  Fig.  1  with  following  equations: 


Ki 


XjAj 


(2) 


Q  =  PtSAe  i  (3) 

where  Ki  and  C,  are  the  thermal  conductance  and  thermal  capacity 
associated  with  each  node  i. 

The  governing  equations  for  p-  and  n-type  thermoelectric 
elements  involve  three  basic  effects,  Seebeck  effect,  Peltier  effect 
and  Thomson  effect.  Besides,  there  are  two  accessorial  effects,  Joule 
effect  and  Fourier  effect.  Gou  [12]  explained  that  Peltier  heat  gen¬ 
erates  on  the  two  sides  of  semiconductor  and  Joule  heat  can  be 
regarded  as  flowing  equally  to  the  two  sides  of  conductors,  so  they 
can  be  disposed  as  a  boundary  heat  fluxes.  Thomson  heat  is  small 
enough  to  be  neglected  since  it  is  a  second-round  effect.  So  the  TEG 
inner  unsteady-state  heat  transfer  equation  can  be  translated  also 
as  Eq.  (1). 
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The  boundary  condition  applied  at  the  hot  plate  surface,  y,  =  0, 
corresponds  to  the  applied  heat  from  electrical  heater  Q,n.  On  the 
top  surface,  y10  =  ech,  the  convection  boundary  equation  is 
employed.  The  boundary  conditions  applied  to  solve  the 
unsteady-state  heat  transfer  Eq.  (1)  for  the  system  shown  in 
Fig.  1  are  listed  below. 

Between  heater  and  hot  plate  surface 


Q-in  =  -•W 


OT| 


-vwfl 


+  'r  T  -  “/ltcl-c',tcl'cay|,6.0 

Between  thermal  compound  and  ceramic  plate 

L 

Between  ceramic  plate  and  thermal  compound 


I2Rpn  =  -/ltd  .cAwiM 


(4) 


Between  hot  plate  and  thermal  compound 

,  .  ari  ,  .  oti 

4Apayjy]^  - 

Between  thermal  compound  and  ceramic  plate 

=  -*ce,hAce<h—^ 

Between  ceramic  plate  and  thermal  compound 


Between  thermal  compound  and  thermoelectric  element 


Between  thermoelectric  element  and  thermal  compound 


where  ;.bp,  '-cp,  Xtci,  he,  hc2,  hb  and  Xpn  are  thermal  conductivities  of 
hot  plate,  cold  plate,  external  thermal  compound,  ceramic  plate, 
internal  thermal  compound,  water-cooling  block  and  thermoelec¬ 
tric  generator,  respectively;  ehp,  ecp,  etci,  ece,  etC2,  eCb  and  epn  are 
lengths  of  hot  plate,  cold  plate,  external  thermal  compound,  cera¬ 
mic  plate,  internal  thermal  compound,  water-cooling  block  and 
thermoelectric  generator,  respectively;  Ahp,  Acp,  Atcl ,  Ace,  Atc2,  Acb 
and  Ap „  are  cross  sectional  areas  of  hot  plate,  cold  plate,  external 
thermal  compound,  ceramic  plate,  internal  thermal  compound, 
water-cooling  block  and  thermoelectric  generator,  respectively; 
Tpn,h  and  Tpnc  are  temperatures  of  hot  and  cold  side  of  semiconduc¬ 
tor,  respectively;  Tcb  and  Tcoo  are  cold  side  temperature  of  water¬ 
cooling  block  and  cooling  fluid  temperature,  respectively;  hcb  is 
coefficient  of  convective  heat  transfer;  /  is  electrical  current  flowing 
through  the  thermoelectric  circuit;  Rpn  and  apn  are  the  electric  resis¬ 
tance  and  Seebeck  coefficient  of  TEG,  respectively;  n  is  the  number 
of  semiconductor  thermocouples;  ocn,  y.p,  Xn,  /,p,  An,  Ap,  e„,  ep,  a„,  op, 
are  the  Seebeck  coefficients,  thermal  conductivities,  cross  sectional 
areas,  lengths  and  electrical  resistivities  of  p-  and  n-type  elements. 

Once  heat  transfer  equation  system,  Eqs.  (4)-(14)  and  addition¬ 
ally  Eqs.  (15)— (19),  is  solved  using  implicit  finite  difference 
method,  the  nodal  temperatures  on  hot  and  cold  side  of  the 
semiconductor  are  found  for  every  time  step.  Consequently,  the 
electrical  and  thermal  outputs  can  be  obtained  every  time  step. 
According  to  Seebeck  effect  the  open  circuit  electromotive  force 
generated  is  expressed  as  follows 


Vo  c(t)  =  <Xp„(Wt)  -  W»)) 


(20) 


The  output  current  /  and  output  voltage  VRf  generated  depends 
on  the  load  resistance  RL  and  it  can  be  calculated  following  Eqs. 
(16)  and  (18),  respectively. 

^(t (21) 
Rpn  +  Rl 

Besides,  short  circuit  current  also  can  be  obtained 


Between  thermal  compound  and  cold  plate 


OTI 


,  OTI 


Between  cold  plate  and  water-cooling  block 

= ~XcbAcb^[K=o 

Between  water-cooling  block  and  cooling  fluid 

With 

ocpn  =  n(ap  -  a„) 

Xpn  =  n(Xp  +  Xn) 

Apn  =  n(Ap  +  An) 

opep  <jne„s 


Rpn  =  n^. 


'  A„ 


_  <Xpn(Tpnlh  -  Tpn,c) 


So  the  output  power  when  RL  is  connected  is  given  by 

P3) 


Deriving  the  Eq.  (23)  with  respect  to  load  resistance,  PR1/J?1  =  0,  it 
can  be  observed  that  maximum  output  power  will  be  reached  when 
load  electric  resistance  is  equal  to  the  internal  electric  resistance. 
Therefore,  substituting  R/.  =  Rpn  into  Eq.  (23)  the  maximal  output 
power  is 

_  0 2pn(Tpn,h(t)  -  Tp„c(t))2 


Pma*(t)=- 


4RP„ 


(24) 


(15) 

(16) 

(17) 

(18) 


The  heat  absorption  on  hot  side  and  heat  release  on  cold  side  of 
semiconductor  are  given  by  Qj,  and  Ci,  respectively. 

Qh(t)  =  Kpn(Tpn,h(t)  ~  WM)  +  MW  WM  “  (25) 


Qc(t)  =  Rpn(Tpnih(t)  -  Tp„c(t))  +  apnI(t)Tp„c(t)  +~I2(t)Rpn 


(26) 


Finally,  the  instantaneous  and  maximal  system  efficiency,  >;  and 
max,  are  calculated  by  Eqs.  (27)  and  (28),  respectively. 


(Rpn  +RL)2Km(Tpn*(t)  -  Tpn,c(f))  +  (Rpn  +Rl)20Cp„I  (1)^(1) -^Rpn +RLfl2  (t)Rpn 

(27) 


(19) 
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= _ ^(WQ-Wt))2 _ 

4Rp„Kpn(Tpn  h(t)  -  Wt))  +4RpnocpnI(t)Tpnjl(t)  -  2RpnI2(t)Rpn 
(28) 

2.3.  Implementation  in  TRNSYS  environment 

The  TRNSYS  software  package  has  been  used  extensively  for 
thermal  system  analysis.  It  has  a  modular  structure  and  consists 
of  individual  subroutines  that  represent  real  physical  devices  or 
utility  components.  The  components  can  be  connected  together 
to  form  complex  systems. 

TRNSYS  uses  a  modular  approach  to  solve  thermal  energy 
systems.  Basically,  it  is  an  equation-solving  program  based  on 
standard  numerical  techniques,  which  provides  different  subrou¬ 
tines  to  find  analytical  solution  to  the  set  of  equations.  However, 
only  a  first-order  linear  differential  equation  solver  is  provided. 
Therefore,  a  subroutine  for  solving  nonlinear  differential  equations 
of  Section  2.2  using  the  finite  difference  and  the  Newton-Raphson 
methods  is  programmed  in  Fortran. 

Fig.  2  shows  the  TRNSYS  project  created  to  simulate  the  TEG 
model.  To  do  so  the  following  input  connections  has  to  be 
established: 

•  The  applied  heating  rate  Qi„,  obtained  from  heater 
specifications. 

•  The  cooling  fluid  temperature,  extracted  from  experimental 
data  measurements. 

•  The  load  resistance  RL,  which  is  used  to  calculate  the  output 
data  when  specific  resistance  is  connected  between  positive 
and  negative  connections.  It  may  be  assigned  or  leave  as  default 
if  no  load  is  considered. 

•  Simulation  time  step  size,  which  is  selected  based  on  the 
required  accuracy  of  the  output  data.  It  also  affects  to  the 
required  computational  time. 

Additionally,  the  geometric  data  and  the  material  properties  of 
the  elements  of  the  thermoelectric  system  studied  have  to  be 
included  into  the  TEG  model  as  TRNSYS  parameter  data.  All  these 
parameters  can  be  shown  in  Tables  1  and  2. 


After  the  simulation,  the  output  values  obtained  from 
Eqs.  (20)— (28)  are  shown  through  built-in  online  plotter. 

3.  Experimental  setup 

To  compare  simulated  results  with  experimental  data,  the 
simulation  must  have  the  same  operating  conditions  as  the 
experiment.  The  validation  of  the  TRNSYS  Type  developed  in  this 
study  uses  experimental  data  from  an  especially  designed  equip¬ 
ment  to  perform  this  test.  The  experimental  setup  allows  establish¬ 
ing  parameters  required  for  the  simulation  of  the  performance  of 
the  TEG,  such  thermal  and  electrical  transients  of  different 
magnitudes. 

Fig.  3  shows  the  schematic  of  the  experiment  with  a  TEG  on  top 
of  a  heat  source.  The  source  temperatures,  T,  and  Th,  are  measured 
by  two  K-type  thermocouples  with  diameter  of  2  mm  placed,  as 
shown  in  Fig.  1,  underneath  the  TEG  and  ceramic  substrate,  respec¬ 
tively.  The  thermocouple  is  sealed  by  conductive  thermal  com¬ 
pound  to  assure  uniform  heat  distribution  on  the  TEG  hot 
surface.  A  water-cooled  block  cools  the  top  surface  of  the  TEG 
through  a  pumped  liquid  cooling  loop  consisting  of  a  pump,  a  heat 
exchanger  and  liquid  lines.  The  pump  circulates  the  fluid  in  the 
loop,  which  picks  up  the  heat  in  the  cold  plate  and  dissipates  it 
through  the  heat  exchanger.  Two  ceramic  plates  made  of  alumina 
(A1203)  are  sandwiching  the  thermo-elements.  The  temperatures  of 
the  top  surface  (cold  side)  of  the  TEG,  T2  and  Tc,  are  also  measured 
by  two  K-type  thermocouples  of  2  mm  mounted  above  ceramic 
substrate  and  cold  plate.  The  hot  and  cold  plates  are  treated  as  heat 
spreaders  to  provide  a  uniform  thermal  field  to  both  hot  and  cold 
sides  of  TEG  module. 

Initially,  the  system  is  at  room  temperature,  measured  as 
16.17  °C,  and  the  electrical  heater  is  switched  on  at  its  first  heating 
level  of  Qin  =  23.6W  (Case  1).  During  temperature  increase,  an 
external  switch  intermittently  connects  and  disconnects  an  electri¬ 
cal  load  of  2.5  Cl  (every  10  s).  This  allows  us  to  obtain  loaded  and 
open-circuit  outputs  simultaneously.  The  transient  temperatures 
at  the  hot  and  cold  surfaces  of  the  TEG,  the  output  voltage  and 
the  output  current  are  recorded  simultaneously  into  a  digital  com¬ 
puter  using  a  data  acquisition  system  (National  Instruments-cRIO). 
The  data  acquisition  is  set  at  a  rate  of  1  s  for  about  3500  s  at  which 
time  the  system  attains  the  steady  state.  Once  system  reaches  the 


Fig.  2.  TRNSYS  project  created  for  TEG  system  I 
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Thermoelectric  parameters 


Parameters  Case  1  Case  2  Case  3  Unit  Source 


23.60  33.50  53  W 

16.45  17.34  18.20  °C 

2.5  2.5  2.5  Q 

163.546  167.273  171.659  pV/K 

-178.098  -188.726  -201.897  pV/K 

1.510  1.453  1.415  W/mK 

1.954  1.925  1.792  W/mK 

10.896  12.001  13.767  pQm 

11.196  13.100  15.910  pQm 


Manufacturer 

Manufacturer 

Manufacturer 

Manufacturer 

Manufacturer 

Manufacturer 


In  Table  1,  the  Seebeck  coefficient,  thermal  conductivity  and 
electrical  resistivity  of  the  p-  and  n-type  elements  are  listed  at 
steady-state  average  temperatures  of  T  =  (Th  +  Tc)/ 2,  which  are 
extracted  from  experimental  data.  All  these  three  parameters  are 
obtained  from  manufacturer’s  datasheet. 

Experiments  are  conducted  under  following  conditions:  the 
temperature  of  cooling  fluid  stays  between  16.45  °C  and  18.20  °C, 
the  air  temperature  is  16.17  °C  and  the  flow  rate  of  the  cooling  fluid 
is  72  1/h. 


4.  Results  and  discussion 


Table  2 

Properties  of  materials  used  in  TEG  tests. 


Source 


7700 

7700 

2700 

2040 

3975 

2040 

2700 


W/mK 
W/mK 
W/mK 
W/mK 
W/mK 
W/mK 
kg/m3 
kg/m3 
kg/m3 
kg/m3 
kg/m3 
kg/m3 
kg/m3 
kg/m3 
J/kgK 
J/kgK 
J/kgK 
J/kgK 
J/kgK 
J/kgK 
J/kgK 
J/kgK 
W/m* 2  K 


Measured 

Measured 

Measured 

Measured 

Measured 

Measured 

Measured 

Estimated 

Manufacturer 

[13] 

Manufacturer 


Manufacturer 


Manufacturer 

Estimated 

Estimated 

Manufacturer 


steady  state,  the  output  voltage  and  current  are  recorded  simulta¬ 
neously  while  load  resistance  varied  from  Ofl  to  100 Cl.  Table  1 
shows  the  initial  values  of  the  test  conducted.  The  experiment  is 
repeated  for  two  additional  hot  side  heat  fluxes  Qjn  of  33.5  W  (Case 

2)  and  53  W  (Case  3)  in  order  to  assess  the  model  at  three  different 
temperature  gradients. 

To  assess  the  TEG  behaviour,  several  thermoelectric  properties 
should  be  tabulated  as  a  reference.  The  parameters  are  determined 
either  according  to  direct  measurement,  from  manufacturer’s  data 
or  from  the  existing  literature  [12,18],  The  geometric  parameters 
and  properties  of  hot  and  cold  plates  and  water-cooling  block  are 
estimated  considering  that  materials  used  are  aluminium  and  coo¬ 
per,  respectively.  Tables  1  and  2  show  the  properties  of  materials 
used  in  the  experimental  setup. 


Based  on  transient  and  steady  state  experiments  and  TRNSYS 
simulations,  the  results  and  a  discussion  about  their  comparison 
is  presented  in  this  section. 


4.1.  Effects  of  step  changes  of  supplied  heat  flow  rates 

The  dynamic  response  characteristics  of  temperature  difference 
between  hot  and  cold  sides  AT  =  Th-Tc,  hot  side  inlet  heat  flux  Qh, 
electric  output  power  Pm  and  closed-circuit  voltage  VKI  with  step 
changes  of  supplied  heat  flow  rate  (&„  are  described  in  Fig.  4.  Both 
the  transitory  experimental  and  simulation  results  are  plotted  ver¬ 
sus  the  transient  time  in  the  same  plot.  The  conditions  imposed  in 
the  experiment  are  listed  in  Table  1. 

During  the  transient,  the  temperature  difference  across  the 
TEG  device  increased  until  the  steady  state  is  achieved.  It  can 
be  observed  that  the  electric  output  characteristic  is  mostly 
dependent  on  the  temperature  difference.  The  hot  side  tempera¬ 
ture  Th  achieved  is  99  °C,  140  °C  and  216  °C  for  cases  1,  2  and  3, 
respectively.  Accordingly  to  experimental  data  obtained  and  Eqs. 
(20)-(28),  the  higher  the  temperature  difference  across  the  TEG 
is,  the  more  electrical  power  generated  and  the  greater  output  sys¬ 
tem  performance  becomes. 


4.2.  Effects  of  step  changes  of  load  current  at  constant  supplied  heat 
flux 

In  Fig.  5,  the  step  changes  effects  of  load  current  under  condi¬ 
tions  of  Table  1  can  be  observed.  The  system  has  been  brought  to 
steady-state  temperatures  of  Fig.  4(a).  It  is  very  interesting  to  note 
the  effect  that  the  load  current  has  on  the  temperature  gradient 
across  the  device,  which  is  the  same  obtained  in  Ref.  [15],  The  load 
currents  imposed  for  cases  1-3  are  0.49  A,  0.71  A  and  1.01  A, 
respectively. 

Every  increase  in  I  leads  to  a  decrease  in  temperature  difference 
AT  due  to  the  greater  value  of  the  Peltier  term  of  Eqs.  (8)  and  (9). 
This  current  increase  also  produces  an  increment  of  hot  and  cold 
side  heat  fluxes.  However,  as  explained  in  Section  2.2,  the  maxi¬ 
mum  output  power  and  efficiency  happens  only  when  load  resis¬ 
tance  is  matched  with  internal  resistance. 

As  it  can  be  appreciated  from  Fig.  5,  the  simulation  results  show 
good  agreement  with  the  experimental  results,  accurately  tracking 
the  electro-thermal  coupled  effects  occurring  in  the  TEG  system. 
The  plotted  curves  of  Th  and  Tc  show  that  steady-state  are  not 
achieved  with  step  changes  of  10  s.  It  also  can  be  shown  that 
parameter  variations  with  load  current  increase  with  higher  values 
of  supplied  heat  flux.  The  low  accuracy  observed  in  Fig.  5(c)  and  (d) 
may  be  caused  by  convection  heat  losses  on  hot  and  cold  plates. 


4.3.  Effects  of  load  resistance  variation  at  steady-state 

The  variation  of  performance  at  steady-state  with  load  resis¬ 
tance  is  shown  in  Fig.  6.  Three  different  conditions  corresponding 
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Fig.  3.  Experimental  test  scheme. 


to  described  cases  in  Table  1  are  selected.  The  system  has  been 
brought  to  the  same  steady-state  temperatures  of  Fig.  4(a). 

Fig.  6  demonstrates  that  the  electrical  output  is  strongly  cou¬ 
pled  to  both  the  thermal  behaviour  and  the  load  resistor  value. 


The  electric  load  resistance  that  produces  the  maximum  electrical 
power  varies  with  the  supplied  heat  flux.  The  power  generation 
and  system  efficiency  are  the  highest  when  the  load  resistance 
value  is  close  to  the  internal  resistance  RL  =  Rpn,  which  is  measured 
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Fig.  5.  Correlation  of  experimental  and  simulated  behaviours  of  the  TEG  system  under  step  changes  of  load  current,  (a)  Hot  side  temperature,  (b)  Cold  side  temperature,  (c) 
Hot  side  heat  flux,  (d)  Cold  side  heat  flux,  (e)  Closed-circuit  voltage,  (f)  Electric  output  power. 


to  be  about  2.5  fi.  As  can  be  shown  in  Fig.  6(d),  the  generated 
power  increases  to  a  maximum  but  then  decreases  to  reach  a 
steady  state.  The  electrical  resistances  at  maximum  power  points 
are  2.4  Cl  (0.65  W),  2.45  n  (1.28  W)  and  2.9  fi  (2.6  W)  for  cases 


1-3,  respectively.  Those  effects  can  be  observed  in  the  theoretical 
model  as  well. 

Besides,  in  Fig.  6(f)-(h)  it  can  be  observed  that  both  the  hot  and 
cold  sides  heat  fluxes,  Qi,  and  Qc,  decrease  with  resistance  until 
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steady  state  is  achieved.  It  demonstrates  that  when  load  resistance 
RL  increase,  the  load  current  I  falls  to  minimum,  which  reduces  the 
Joule  and  Peltier  effect  of  the  Eqs.  (25)  and  (26).  Consequently, 
when  load  current  is  very  low,  those  effects  are  despicable  and  it 
is  only  the  Fourier’s  law  that  governs  the  heat  transfer  equation. 

4.4.  Model  uncertainty 

The  TEG  model  developed  in  this  work  is  validated  by  compar¬ 
ing  the  experimental  data  with  the  numerical  data  predicted  by 
our  model  at  the  same  conditions.  The  calculus  of  the  error  is  per¬ 
formed  by  the  comparison  of  the  temperature  difference  and  the 
electrical  output  power  of  each  point  predicted  by  the  model  with 
the  same  point  measured  experimentally.  The  value  of  RMSE  is 
4.85  °C  and  0.0648  W  for  temperature  difference  and  electrical 
output  power,  respectively.  The  model  presents  a  normalized  root 
mean  square  error  of  3.53%  and  2.33%  for  temperature  difference 
and  electrical  output  power,  respectively. 

Although  the  model  presented  has  slightly  higher  temperature 
difference  RMSE  and  NRMSE  than  the  model  developed  by  Monte- 
cucco  [15]  (3.58  °C  and  2.76%,  respectively),  the  output  power  pre¬ 
diction  is  significantly  better  (0.2  W  and  4.55%,  respectively).  It  is 
important  to  note  that  the  model  accuracy  is  markedly  affected 
by  the  thermoelectric  parameters  introduced  into  the  simulation 
model.  Consequently,  the  more  exact  the  thermoelectric  parame¬ 
ters  are,  the  more  precise  the  outputs  will  be. 

5.  Conclusions 

Although  thermoelectric  phenomena  have  been  used  for  heat¬ 
ing  and  cooling  applications  quite  extensively,  electricity  genera¬ 
tion  has  only  seen  very  limited  market  in  niche  applications  and 
it  is  only  in  recent  years  that  interest  has  increased  regarding 
new  applications  of  energy  generation  through  thermoelectric  har¬ 
vesting.  The  widespread  use  of  TEGs  depends  on  its  optimization. 

The  TRNSYS  component  presented  in  this  paper  has  been  devel¬ 
oped  for  this  purpose  and  is  described  and  validated  using  experi¬ 
mental  data.  The  proposed  component  is  able  to  cope  with  thermal 
and  electrical  dynamics.  The  comparison  of  results  between  theo¬ 
retic  analysis  and  experiment  has  approved  the  reasonability  of 
the  new  component. 

On  the  other  hand,  the  TRNSYS  simulation  runs  without  inter¬ 
ruptions  or  delays,  therefore  the  numerical  model,  and  conse¬ 
quently  the  new  component,  can  be  considered  well  optimized. 
Therefore,  this  system  model  can  be  used  in  performance  optimi¬ 
zation  and  further  application  of  thermoelectric  generation. 
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